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PREFACE 


The  Project  RAND  research  program  consists  in  part 
of  basic  supporting  studies  in  rnafnematics.  One  aspect 
of  this  involves  the  investigat±on  of  the  approximate 
behavior  of  solutions  of  differential  equations. 

In  the  present  Memorandum  the  author  discusses  the 
problem  of  determining  useful  analytic  approximations  bo 
solutions  of  equations  of  the  form 

L(u)  +  (a(p)  +  >b(p))T_i  =  0 
vdiere  L(u)  is  a  linear  differential  operator. 
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SUMMARY 


A  problem  of  continuing  interest  is  that  of 
obtaining  approximate  solutions  of  the  functional 
equation 

L(u)  +  (a(p)  -f  Xb(p))u  =  0, 

v;here  L  is  a  linear  transformation,  in  terms  of  the 
solution  of  the  unperturbed  equation 

L(u)  -1-  a(p)u  =  0. 

Using  the  Green's  function,  or  equivalent  techniques, 
and  regarding  the  term  involving  X  as  a  forcing  term, 
we  can  convert  the  first  equation  to  the  form 

u  =  f  -t-  XT(u) , 

where  T  is  a  linear  transformation. 

We  present  a  new  approach  to  problems  of  this 
nature  using  the  classical  technique  of  continued 


fractions. 
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A  NH-J  EORl'IALISM  IN  PERTURBATION  THEORY 
lUSING  CONTINUED  FRACTIONS 


1 .  INTRODUCTIQ!^ 

A  problem  olC  continuing  Interest  is  that  of  obtain¬ 
ing  approximate  ssolutions  of  the  functional  equation 

(1.1)  L(u)  4  (a(p)  4-  Xb(p))u  =  0, 

I'/here  L  is  a  linear  differential  operator,  in  terms  of 
the  solution  of  fthe  unperturbed  equation 

(1.2)  L(u)  4  a(p)u  =  0. 

Using  the  Green 'ss  function  of  (1.2),  or  equivalent 
techniques,  and  iregarding  the  term  involving  X  as  a 
forcing  term,  w€  can  convert  (1.1)  into  an  equation  of 
the  form 

(1.3)  u  =  f  XT(u), 

•where  T  is  a  linear  transformation. 

We  shall  present  a  new  approach  to  problems  of  this 
nature  using  the;  classical  technique  of  continued 
fractions.  In  subsequent  papers  we  shall  discuss  the 
rigorous  aspectsi,  and  the  value  of  this  approach  in 
obtaining  rational  approximations  outside  the  circle  of 
convergence  of  the  Li ouvi lie— Neumann  solution  of  (1.3) • 
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2.  CONTINUED  FRACTION  EXPANSION 
From  (1.3)  we  have 

(2.1)  T^(u)  =  T^(f)  +  XT^'^^(u),  n  =  0,1,2, - 

Hence,  combining  the  n— th  and  (n  +  l)— st  expressions. 


a  three— tenn  recurrence  relation.  Thus 


(2.3) 


XT^(f  )  ) 

T^+^(u)/T^+^(u)  ' 


Thus,  formally  (setting  T^(f)  =  Tj^)^ 


(2.4) 


u  _  "^'o  ^  “N  ^'^c/’^i 

TTTTT  =  (^  TT  +  ^  j  -  VTT-^  XT  /T 

-  TTTa 


Retuming  to  (1.3),  we  have  T(u)  =  (u  —  f)/X.  Combining 
this  result  with  (2.4),  we  have  the  follox^ng  representa¬ 
tion  for  the  solution: 


(2.5)  u  = 


''VA 


XToTg 

TTJ+XTg) 
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4 .  EXAMPLES — II 

Consider  the  case  where  we  replace  X  by  —  X, 
f  >  0  and  T  is  a  positivity-preserving  transformation. 
Then  a  fonnal  solution  of  u  +  XT(u)  =  f  has  the  form 


(4.1)  u  = 


Hence  If  -  X  >  0  Tor  all  n,  we  see  that  all 

of  the  convergents  to  (4.1)  are  obviously  nonnegative. 
We  thus  obtain  a  simple  sequence  of  lower  and  upper 
bounds  to  the  solution  of  u  +  XT(u)  =  f. 


5 .  EXAMPLES - III 

Interesting  analogues  of  asynptotic  expansions  are 
obtained  upon  considering  equations  of  the  form 

(5.1)  u  =  f  +  e  ^  , 

dt 

and  the  matrix  analogues. 


